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cond-niat/0610469 , to calculate the anisotropy effect in the interaction of two level 
systems with phonons in disordered crystals. We particularize our calculations to 
cubic crystals and compare them with the available experimental data to extract 



_i4 ^ the parameters of the model. With these parameters we calculate the interaction 

Q ' of the dynamical defects in the disordered crystal with phonons (or sound waves) 



propagating along other crystalographic directions, providing in this way a method 
to investigate if the anisotropy comes from the two-level systems being preferably 
oriented in a certain direction or solely from the lattice anisotropy with the two-level 
systems being isotropically oriented. 

PACS numbers: 61.43.-j,62.40.+i,63.20.kp 

INTRODUCTION 

The low temperature acoustic and thermal properties of amorphous, glassy materials 
are remarkably similar and they can be explained to a large extent by assuming that the 
material contains a large number of dynamic defects. These dynamic defects are tunneling 
systems (TS) and are modeled by an ensemble of two- level systems (TLS) [l|, |2]. Crystals 
with defects-with a large enough amount of disorder-exhibit also glass-like properties, but 
these properties are not so universal and, even more, they are not isotropic-like it is, for 
example, the sound absorbtion and velocity change, which depends on the crystalographic 
direction in which the sound propagates 
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Since a detailed microscopic model of tunneling systems in glassy materials is still not 
available, the study of disordered crystals is especially interesting because it offers an addi- 
tional opportunity for their clarification: in some materials we know quite well which are the 
entities that tunnel between different equilibrium positions. Beside this, the anisotropy of 
the TLS-sound wave interaction in crystals represents another challange to the interaction 
models of TLSs which requires clarification. 

In this paper we give an explanation for the anisotropy observed in the glass-like properties 
of general, disordered crystals, by employing a model recently published ^. In this model 
we assume that each TLS is characterized by a direction in space, call it t-this might be the 
direction defined by the two potential wells of the tunneling system, or the axis of rotation 
of the tunneling entity-and we introduce a coupling between the TLS and a strain field [S], 
which is dependent on the amplitude of [S] at the position of the TLS and on the orientation 
of t with respect to [S] . In (^ the model was applied to an amorphous solid, assuming that 
the directions t are isotropically distributed, and the effective coupling of an elastic wave 
with a TLS was calculated as the average over the directions of the TLS. In this way it was 
proven on very general grounds that, on average, the longitudinal waves couple with the 
TLSs stronger than the transversal waves-in standard notations, 7« > 7t |4|. 

In a disordered crystal with TLSs, there could be at least two sources of anisotropy. The 
first one is that the TLSs might not be anymore isotropically oriented so the effective coupling 
of elastic waves with them depends, through t, on the waves direction of propagation and on 
their polarization. The second source of anisotropy is that besides the relative orientation of 
t and [S] , the symmetry of the crystal is manifested also in the interaction of elastic waves 



with TLSs 
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This leads to anisotropy effects in the interaction of elastic waves with 
TLS even if the TLS distribution is isotropic. 

In this paper we shall analyse mainly the second type of anisotropy. We shall assume 
that the TLS orientations are isotropically distributed and we shall calculate the anisotropy 
effects imposed only by the lattice symmetries onto the interaction hamiltonian. Existing 
experimental data (see Ref. |6| and references therein) will enable us to obtain relations 
between the parameters of the model for a cubic lattice. Using these parameters we can make 
predictions about the absorbtion of elastic waves propagating along other crystalographic 
directions or having different polarizations. These predictions could constitute a first test 
for the isotropy of the TLS orientations in a specific crystal. 



ANISOTROPIC INTERACTION OF TWO-LEVEL SYSTEMS WITH SOUND 

WAVES 

In the standard tunneling model (STM), the Hamiltonian of an isolated TLS is written 



in a two-dimensional basis as 
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M A A _ 1 A -A 

i^TLs = ya. - -a. = 2 1^ _^ _^ I (1) 

where A is called the asymmetry of the potential and A the tunnel splitting. The basis in 
which is written the Hamiltonian of Eq. ([1]) is chosen in such a way that a perturbation to 
the TLS, caused by a strain field, say [S], is described by a diagonal Hamiltonian, 
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with (5 = 2(7] : [S] and [7] a second rank tensor of coupling constants; by ":" we denote the 
dyadic product. Typically, in the STM one considers the coupling of TLSs with transversally 
or longitudinally polarized sound waves, so not too much attention has been given to the 
[7] tensor and in general 6 is written simply as 6 = 2'~fi^tSi^t, with 7 and S being scalars {S is 
the amplitude of the strain field) and / and t denoting the longitudinal (/) or the transversal 
(t) polarization of the sound wave, respectively. Such a simple description of the TLS-strain 
field interaction has several short comings-e.g. 6 is not invariant and even leads to physical 
ambiguities at the rotation of the coordinates axes-and cannot account for the different 
couplingof the TLS with sound waves propagating in different directions. In consequence, 
in Ref. (j] it is proposed a model which eliminates the shortcomings and takes into account 
the symmetries of the material in which the TLSs are embeded and the orientation of the 
TLS with respect to the strain field. Let us describe briefly how this is done. 

We construct form the components oft the simple 3x3 symmetric tensor, [T], of compo- 
nents Tij = titj and we introduce the forth rank tensor of TLS-strain field coupling constants, 
[[/?]]. With these two objects, we build the general tensor [7], as •jij = Tfci_Rfc/jj-throughout 
this paper we assume summation over the repeated indices. The forth rank tensor [[/?]] 
has a similar structure as the forth rank tensor JTcjl of stifness constants and reflects the 
symmetries of the crystal that contains the TLS 4, LZl]. 

For the convenience of the calculations we work here, like in Refs. y, la, LZ|, in ab- 
breviated subscript notations and write [T] and [S] as the six elements vectors, T = 



(Tn,T22,T33,2T23,2Ti3,2Ti2)* and S = (^n, ^22, ^33, 2^23, 2^13, 2^12)*, where by "■*" we 
denote the transpose. Following the notations of Auld |8|, the components of the symmetric 
tensors will be denoted in abbreviated subscript notations by a single, upper case subscript- 
e.g. Tj, Si, and T3 = T33 = t|; also in abbreviated subscript notations, the tensors [[R]] and 
[[c]] will be written as 6 x 6 matrices, [R] and [c], of components i?rj and cjj, respectively. 



Q,y,i7i. 



Putting all these together we get the expression 6 = 2T* ■ [_R] ■ S 

Having now the full expression for the interaction Hamiltonian, Hi, we can calculate 
the amplitude of excitation of a TLS, of parameters A and A, by a phonon of wave-vector 
k and polarization a; we denote by n]f_cr the number of phonons on the mode (k, a) after 
the TLS excitation process. The displacement field of the phonon, Uk^, is normalized to 
A'ko- = Jh/{2VpuJka-), and has the strain field Sko- = V^Uko- (where by V5 we denote the 
symmetric gradient). This way we get 

(nk., T l^il^k. + 1, i)u,. = --V^T' ■ [R] ■ Sk. (3) 

where e = v^A^~+~A^ is the excitation energy of the TLS. Therefore the phonon scattering 
rate by a TLS in the ground state is 
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\T'-[R]-S^^\'6{e-hu). (4) 



The main characteristic of the TLS-elastic strain interaction, is contained in the quantity 
Mk,o-(t) = T* ■ [R] ■ Skcr, which we shall calculate next. 

As mentioned above, [R], like [c], reflects the symmetries of the lattice. The most general 
type of lattice is triclinic, in which case [R] is symmetric and contains 21 independent 
constants. Such a lattice is very complex and in general does not sustain simple transversally 
or longitudinally polarized elastic waves, but instead, the elastic waves propagating through 
the crystal will be complex superpositions of longitudinally and transversally polarized plane- 
waves. So, to start with a simpler case and also to be able to compare our calculations with 
available experimental data 6|], we shall focus in this paper on lattices with cubic symmetry. 
The tensor [R] for the cubic lattice is very similar to the one for an isotropic material 4, l5|, LZl] , 
but it contains 3 independent constants instead of 2, like in the isotropic case. So we can 



preserve the notations of Refs. 4, l5|, l3], and write 
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without imposing the isotropy constraint, ( + 2^ = 1; similarly, the tensor of elastic stiffness 
constants is 
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Using [c] we can write the Christoffel equation to find u and S for the elastic waves propa- 
gating in different directions and then we can calculate M for any t, A and A. In the end, we 
average over the ensemble of TLSs, to determine the attenuation of the elastic wave or the 
scattering rate of the phonon. We shall apply this procedure for strain fields corresponding 
to elastic waves propagating along the crystalographic directions (100), (110) and (111) of 
the cubic lattice. Along these directions, the cubic lattice can sustain simple, longitudinally 
and transversally polarized elastic waves, for any allowed values of the parameters Cn, C12, 
and C44. 

Solving the Christoffel equations we find that the sound velocities of the longitudi- 
nal waves propagating in the (100), (110) and (111) directions are c;^(ioo) = yCn/p, 
Q,(iio> = \/(cii + C12 + C44)/p, and q,(iii) = y^(cii + 2ci2 + 2cm)/p, respectively. Similarly, 
the sound velocities of the transversal waves propagating in the (100) and (111) directions are 
Ct, (100) = Y C44/P and Ct^(iii) = \/(^ii ^ '^12 ^ <^aa)/ Pi respectively, whereas for the transver- 
sal waves propagating in the (110) direction the sound velocity depends on the direction 
of polarization: if the wave is polarized in the (100) direction (and perpendicular on the 
direction of propagation), the sound velocity is q ,-^1^-. = 



C44/P, and if the wave is polarized 



in the (HO) direction (and also perpendicular on the direction of propagation), the sound 
velocity is cj. ,-^/qv = \/(cii — C12 — c^a)/ P- Now we can calculate M for these three directions 
of propagation. 

Since the three directions, x, y, and z are equivalent, let us take the (100) direction 
as the z direction. We also define t by the angles 9 (nutation) and cf) (precession), as 
t = (sin ^ cos 0, sin 6' sin 0, cos 6')*. With these conventions, we get for the longitudinal wave, 
Ufc,/r) = iVze^'=-^ 

Mui^i = ik^Nki^i[C + cos" ^(1 - 0] (7a) 

(to simplify the expressions without reducing the clarity, we shall always drop the exponential 
from the expressions of M and the subscripts of A^; the implicit subscripts of A^ are always 
the same as the ones of M and u) and for the two, reciprocally perpendicular, transversal 
waves, Ufcz,t,a;(r) = iVxe^^^'"" and vikx,t,yij) = iVye'^^-"", 

Mki,t,c, = ik^^N sm{29) cos{^) (7b) 

Mki,t,y = ik^^N sm{29) sm{(P) (7c) 



For the waves propagating in the (111) direction we get the following results. For the 
longitudinal wave, u^^+^^-^/^/r) = A^^^^^exp 



^k^^ . r 



M,(.+^+,)/^_, = iV^{[2sin(2e)(sin0 + cos0) (8a) 



+2 sin(20) sin^(^)]^ + 2C + 1} 



ik^^^-r 



and for the two transversal waves, '^k(x+y+z)/^/3 1 pi (^) ~ -^^i ^-^P 
Ufc(x+y+z)/V3,t,p2(r) = A^P2exp [ik^^^^ ■ r], with polarizations pi = ^^, ps = ^ 



and 

we 



have 



and 



1 n'~^ 

M kjii+y+i) ^ = N^^[2^{cos9 — cos(f)sm9) sm9sm(j) 

+ (cos2^-cosVsin^^)(l-C)] (8b) 



^.(x+y+z)/V3,t,p. = iV^[sin(2^)(sin0-cos0)e 

-sin2ecos(20)(l-C)] (8c) 



respectively. 



For the longitudinal wave, U;,(-^_,_^w^^^(r) = N^^e ~^ "", propagating in the (110) di- 
rection, 

^Mx+y)/V2,/ = iV^[2sin(20)sin2^e + (l + cos2^)C 

+ sin^ Q\ (9a) 

and for the two transversal waves in the same direction, u^(x+y)/v^jp' (r) = 
N-p'^ exp ik^^^ ■ r and u^(x+y)/y^,ty (r) = Np'2exp ik^^^ ■ r , with polarizations p[ = z, 
P2 = "^^ , we have 



and 
respectively. 



^.(x+y)/^,,pi = ^^ V ""^^ ^ ^O<20)(C - 1), (9C) 



Now we can calculate the phonon's scattering rates, by averaging Fko- of Eq. (j4j) over 
the distribution of TLS parameters. A, A, 6 and and taking into account the scattering 
of phonons from and into the mode (k, a). We assume that the parameters A and A are 
independent of the parameters 6 and 6, and their distribution is the standard P(A, A) = 

n Q 

Pq/A, where Pq is a constant [l|, ^. We change the variables A and A into the variables 



e and u = A/e, with the probability distribution P(e, u) = Pq/ \u\/1 — u^] and we assume 
that the fraction of excited TLSs, of energy e, is thermal and corresponds to a temperature 
T: n^^^^^ = (l + e'^/^'s'^)"^. The distribution over 9 and (p, say f{6, (p), is unknown. Plugging 
all these quantities into the standard scattering rate calculation, we get 

PotanhfT^T^) pn 

r^^ = ^T^^^^k. / sin^rf^ 

Zri Jo 

JO 

Potanh(:rT:^) ^ ^ 

= ^A?^^nk.(|Mk.(t)H. (10) 

Now, if we would know f{9,(j)), we could use the expressions ([7]), ([H]), and (^ for M, to 
calculate scattering rates of phonons propagating in the three different directions. 

Since, as we mentioned in the Introduction, we have no microscopic model for f{6, 0), we 
shall assume that f{9, (p) is constant (i.e. TLSs are isotropically oriented) and comparing 
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our calculations with experimental results, we shall obtain relations between the parameters 
C and ^. Using this assumption and Eq. fl7a|) . we get for the longitudinally polarized phonon 
propagating in the (100) direction, 

^^-' = 15 h '""^ [2k^) ' ^''"^ 

while for the transversally polarized waves, both Eqs. ( 17b[) and ( ITcl) give the same result, 

'^^'* 15 h \2kBTj ' ^ ^ 

where we dropped also the obvious subscripts of the population number, n. 
Similarly, in the direction (111) we get 

1 _ 5 + 20C + 20C^ + 16^2 2nPok^N'^n^'^ 

±y± 



rfe(^+l+^^^ 45 n 



and 



T-' 



X tanh ( ^ ) , (12a) 



2[(1 - C)' + 2^=^] 27TPoeN^n^^ 



^^±^,t 45 n, 

xt^h(5^). (12b) 

where again, the two transversally polarized waves, Eqs. fISbp and fl8c|) . give the same result. 
Finally, for the phonons propagating along the (110) direction we obtain the average 
scattering rates. 



T-' 



2 + 6C + 7(2 + 4^2 27iPoeN^n^^ 



'-^^,1 15 h 



V2 



X tanh ( — — — - ) , (13a) 

\2kBTj ' ^ ^ 



for the longitudinal wave. 



for the transversal wave polarized in the p'^ direction, and 
for the transversal wave polarized in the pg direction. 
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In the STM formalism, with 5 = 2'-)^S„-a = l,t-, the transition rates are [4| 



therefore the Eqs. flTT|) - flT3|) give the expressions for the 7^s for the phonons propagating 



in different crystalographic directions. Notice that if we impose the condition 



or isotropy, 



C + 2^ = 1, all equations ( ITT|) -( TT3l) reduce to the isotropic expressions of Ref. |4|. 

Karen Topp, Robert Pohl, and coworkers (see Ref. 6| and references therein) measured 
7 in the crystalographic directions (100) and (111) of the cubic lattice of Ca stabilized 
Zirconium. They obtained a ratio between 7^ in the (HI) direction and 7^ in the (100) 
direction, of about 1.7. Using this result and Eqs. ( lllbl) and ( 112bl) . we obtain a relation 
between ( and ^: 

(^) =18.4. (15) 

Notice that the lower symmetry of the cubic lattice modifies the relation (1 — C)^/C^ = 4, 
satisfied in an isotropic medium, into the relation ( ITSl) . More experimental data would enable 
one to check this relation for other polarizations or propagation directions and eventually 
even to calculate the elements of [R] . If relation flTSl) does not hold true for any propagation 
direction and polarization of the elastic wave, then the distribution of TLS orientations is 
not isotropic. 

CONCLUSIONS 



a 



We applied the formalism introduced in Ref. [J] to describe the interaction of phonon 
modes (or elastic waves) with the ensemble of two-level systems (TLS) in a disordered cubic 
crystal. We showed that the interaction is anisotropic and-in the language of the stan- 
dard tunneling model-the coupling constants 7/ and 7^ depend on the phonon propagation 
direction. We focused our calculations on phonons propagating along the crystalographic 
directions (100), (110) and (111), for which we gave explicit expressions for the coupling 
constants. Using the experimental results of Topp, Pohl, and coworkers (see Ref. |6| and ref- 
erences therein) we compared the 7iS corresponding to the (100) and (111) crystalographic 
directions and from here we obtained a relation between the the two parameters of the 
model, ( and ^, that describe the anisotropy of the interaction. Nevertheless, more experi- 
mental results are needed (at least •jt in one more direction or a 7^) to fully determine these 



10 

parameters and make prediction about the interaction of the TLS system with phonons 
propagating in any direction. Having these predictions, one then could draw conclusions 
about the isotropy of the TLS orientations in the material. 
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